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> , Abstract 

. We study confined solutions of certain evolutionary partial differ- 

O ! ential equations (pde) in 1 + 1 space-time. The pde we study are Lie- 

■ Poisson Hamiltonian systems for quadratic Hamiltonians defined on the 

. dual of the Lie algebra of vector fields on the real line. These systems 

0> I are also Euler-Poincare equations for geodesic motion on the diffeomor- 

phism group in the sense of the Arnold program for ideal fluids, but 
where the kinetic energy metric is different from the norm of the 
velocity. These pde possess a finite-dimensional invariant manifold of 
particle-like (measure-valued) solutions we call "pulsons." We solve the 
c/3 ■ particle dynamics of the two-pulson interaction analytically as a canoni- 

cal Hamiltonian system for geodesic motion with two degrees of freedom 
and a conserved momentum. The result of this two-pulson interaction 



■ for rear-end collisions is elastic scattering with a phase shift, as occurs 



with solitons. In contrast, head-on antisymmetric collisons of pulsons 
tend to form singularities. 

Numerical simulations of these pde show that their evolution by 
geodesic dynamics for confined (or compact) initial conditions in various 
nonintegrable cases possesses the same type of multi-soliton behavior 
(elastic collisons, asymptotic sorting by pulse height) as the correspond- 
ing integrable cases do. We conjecture this behavior occurs because the 
integrable two-pulson interactions dominate the dynamics on the in- 
variant pulson manifold, and this dynamics dominates the pde initial 
value problem for most choices of confined pulses and initial conditions 
of finite extent. 
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1 Introduction 

We study particle-like dynamics on a set of finite-dimensional invariant man- 
ifolds for certain measure-valued solutions, called "pulsons," of the family of 
evolutionary integral partial differential equations given by 



TJit + urrix + 2mUx = , 



lim m 



0, 



1.1^ 



Here subscripts denote partial derivatives, m:MxM^M is a real, measure- 
valued map on 1 -(- 1 space-time with coordinates x, t, and the function u is 
defined by the convolution integral. 



u{x, t) = g{x - y) m{y, t) dy =: g * m . 

J — oo 



(1.2) 
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The integral kernel, or Green's function, g{x), is taken to be even, g{—x) = 
g{x), of confined spatial extent, and such that the quadratic integral quantity 

1 SH 
H=- mg*mdx with — = u, (1-3) 

is positive definite; so that H defines a norm. Physically, m is a momentum 
density, associated to a velocity distribution u = g * m, and H is the kinetic 
energy for the dynamics. Evenness of the Green's function g{x) implies that 
equation ( |1 . ![ ) conserves the kinetic energy H and the total momentum P = 
I^oo ^ '^^1 solutions m of equation (|1.1|) that vanish at spatial infinity. The 
"mass" M = y/mdx is also conserved for such solutions. 

In a certain formal sense, equation ( |1 . 1|) is hyperbolic, as it follows from a 
pair of equations of hydrodynamic type, namely 

it + ui^ = 0, r]t + {r]u)^ = , with u = -g * (?74) . (1.4) 

This pair of equations implies equation ( p..l|) , upon setting m = —Tjir^. Note 
that i is preserved along flow lines of u and the integrated "mass" / rjdx 
is conserved for solutions of equations (|1.4| ) for which u vanishes at spatial 
infinity. Note also that ix satisfies the same equation as r] does. Thus, these 
equations preserve the relation rj = ±ix, provided it holds initially. Hence, 
equation ( |1 . ![ ) preserves sign-definitiveness of m, provided m can be initially 
expressed as m = ±17^, for some function ri{x, 0) on the real line. Moreover, r/^ 
also satisfies equation ( |1.1D , and this equation is the second in a hierarchy of 
equations implied by ( |1.4| ) for the powers of rj. Because of its special geometric 
properties, we shall concentrate our attention here on the family of equations 
( |1.1| ) for various choices of the Green's function, g{x). 



Hamiltonian and geodesic properties. Equation ( p. . 1| ) is expressible in 
Lie-Poisson Hamiltonian form rrit = {m, H} with Hamiltonian H given in 
equation (|1.3|) and Lie-Poisson bracket defined by 



{r,HY = — / m - — , - — dx = — - — [om + mo) - — dx , (1.5) 
\_om om\ dm dm 

where we have integrated by parts and d denotes the operator d/dx. This 
Lie-Poisson bracket is defined on g*, the dual of the Lie algebra g of vector 
fields on the real line that vanish at spatial infinity and possess the Lie bracket 
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operation ad : q x q ^ q defined by adfh = [/, h] = fh^ — hf^, for f,hEQ. 
Tlius, equation may be rewritten for m G 0* as 



6H 

rrit = {m, H} = -{dm + m9)— = -ad}fj/sm^ > (1-6) 



where ad* is the operation on g* dual to the ati-operation on g. By Legendre 
transforming from g* to g using the relation 6H/6m = u E g, equation (|1.6| ) 
may be re-expressed as an Euler-Poincare equation 0, 



for the Lagrangian 



oo 



1 



oo 



L = I mudx — H = - I uQudx, {l.t 



where g{x) is the Green's function for the self-adjoint positive operator Q, i.e., 

Qg{x)=26{x), (1.9) 

with Dirac measure 5{x). Thus, m = 6L/6u = Qu is the momentum density 
conjugate to the velocity u & g, and u satisfies the Euler-Poincare equation, 

d 

—Qu = —ad*^Qu. (1-10) 

Replacing m by Qu in (|1.1|) reveals that this equation is not Galilean invariant 
(i.e., it changes form under x — x + ct, m — m + c, t — t) unless m — m under 
such a transformation. Typically, m m + k, with k, a constant depending 
on c under a Galilean transformation, so that equation (|1.1| ) becomes 

rrit + uiTLx + 2mUx + 2kUx = , (l-H) 

which for k 7^ introduces linear dispersion. Thus, equation ( |1 . 1|) should best 
be considered as a family of equations parameterized by the function g{x) and 
the constant k, here taken as k = 0. The aim of this paper is to explore 
the solution behavior of this family of equations for the initial value problem, 
under variations of g{x) and the initial conditions. The effects of k 7^ will 
also be discussed briefly in the Appendix. 
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Equation ( p,.10| ) is formally the equation for geodesic motion on the dif- 
feomorphism group with respect to the metric given by the Lagrangian L in 
equation ( p..8| ), which is right-invariant under the action of the diffeomorphism 
group. See Holm, Marsden and Ratiu for detailed discussions, applications 
and references to Euler-Poincare equations of this type for ideal fluids and plas- 
mas. See Ovsienko and Khesin and Segal for discussions of the Korteweg- 
de Vries (KdV) equation from the viewpoint of the Euler-Poincare theory of 
geodesic motion on the Bott-Virosoro group. Misiolek has given a similar 
interpretation to equation (|1.1| ) in the integrable case, g{x) = e"'^' [§], [^. See 
also Holm et al. [Q, |T0| for generalizations of equation ( |1.1| ) to higher spa- 



tial dimensions, and Shkoller [|rT|] and Holm et al. |12] for generalizations to 



Riemannian manifolds. Some of the most interesting solutions of the systems 
we study actually leave the diffeomorphism group due to a loss of regularity. 
Such solutions must be interpreted in the sense of generalized flows, as done 



by Brenier and Shnirelman |[T^. The functional-analytic study of these 
solutions is made in Holm et al. [15|. In this paper we shall formally consider 
such solutions to represent geodesic motion and hence keep the above nomen- 
clature. These dynamics can be formulated either in a periodic domain, or 
on the real line. For the analysis here we shall work on the real line. The 
numerics will be conducted in a periodic domain. 



Completely integrable cases. The Hamiltonian system ( |1 . 1| ) is known to 
be completely integrable in the Liouville-Arnold sense for three cases: g = 
5{x); g = 1 — \x\ on \x\ < 1 (and g = elsewhere); and g = e~'^'. These cases 
are identified with the following Lie-Poisson equations and solution behavior. 

g = 6{x) , rrit + Smrrix = shocks, 

f for X X 

9=\q for|x|>l' i'^t + uu^)^^ = ^{uD^ compactons, 

(7 = e~l^'l , {1 — d'^){ut + uux) = — {u"^ + ^ul)x peakons. 

(1.12) 

The first case yields the Riemann equation, which governs the formation of 
shocks. The second case yields a Galilean invariant equation in the Harry Dym 
hierarchy (at the KdV shallow water position) which describes the propagation 
of weakly nonlinear orientation waves in a massive nematic liquid crystal direc- 
tor field ||16[. The compacton solutions of this equation are triangular pulses 
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with compact support that interact as sohtons. The third case, g = e~'^', 
yields the Camassa-Holm equation, whose peakon solutions describe a lim- 
iting situation of unidirectional shallow water dynamics [Q, 0. Integrable 
compactons and peakons are discussed further in Q, [1^ for motion on 



the real line, and in |]T^, [jTsl and references therein for the periodic case. 



See also references [|T^, [^] for other recent discussions of the Camassa-Holm 
equation. Note that, from ( |1.12| ), the Dym shallow water equation for com- 
pactons arises as a "high wavenumber limit" of the Camassa-Holm equation 
for peakons. The latter two integrable cases each have an associated isospec- 
tral problem whose discrete spectrum determines the asymptotic speeds of 
the solitons that emerge from their initial conditions. Here we shall use these 
integrable solutions as comparisons in assessing the nonintegrable, but quali- 
tatively identical pulse interaction behavior for other choices of g{x), which 
turns out to be the pulse profile. 



Traveling waves. Traveling wave solutions of equation ( |1 . 1| ) of the form 
m{x — ct), with u{x — ct) = g * m, satisfy 

— {c — u)m' + 2mu' = 0, (1-13) 

whose first integral is 

(c — u)^m = = const. (1-14) 

These traveling waves are critical points of the sum of conserved quantities, 
H — cP + 2aM. The boundary condition lim|2;|^oo 1^ = implies that the 
constant a must vanish. Therefore, m must vanish except where g*m = u = c. 
Thus, equation ( |1 . 1|) admits Dirac measure valued (particle-like) traveling wave 
solutions m = 2c6{x — ct), for which u = cg{x — ct). Substitution of these 
formulae for m and u into (|1 . 1|) and integration by parts against a smooth test 
function (using g'{0) = from evenness of g) verifies this solution, provided 
g{0) = 1. Therefore, the traveling wave's speed c is equal to the peak height of 
its velocity profile, and the Green's function g{x — ct) is the normalized shape 
of this profile. We shall study the interaction dynamics of a superposition of 
"pulsons" in the following form 



N 

m{x, ^) = ^ 2pi{t)5{x - qi{t) ) , 

i=l 



(1.15) 
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for which the velocity superposes as [21 

N 



u{x,t) = Y,Mt)9{x-qi{t)), (1.16) 
1=1 

and the Hamiltonian H in ( |1.3|) becomes (up to an unimportant factor of 2) 

1 ^ 



Substitution of these superpositions of particle-hke solutions into ( |1.1|) shows 
that they form an invariant manifold under the Lie-Poisson dynamics of ( |1.1| ), 
provided the parameters {pi,qi), i = 1,...,N, satisfy the finite-dimensional 
particle dynamics equations, 



N N 



Pi = -Pi^Pj9'{qi-qj) , qi = ^Pj9{qi-qj) ■ (i-i8) 



These are precisely Hamilton's canonical equations for the "collective" Hamil- 
tonian Hfyf in ( [LIT] ). These equations describe geodesic motion on an A^- 
dimensional surface with coordinates qi, i = 1, . . . ,N, and co- metric = 
9{qi—qj)- Equations ( |1.18| ) are integrable for any finite N for two cases P3| ; 
namely, g{x) = A -|- /icos(z/x) + /x' sin(z/|x|) and g{x) = A -|- a\x\ + where 
{X, fi, fi',h',a, (3} are free constants. These cases correspond to the pde for 
compactons and peakons in ( |1.12| ) known to be integrable by the isospec- 
tral method, when g{x) is required to be spatially confined. We shall choose 
Green's functions g that are spatially confined (or have compact support); so 
that g{x) and g'{x) are negligible for |a;| > D, with D the interaction range 
between particles. Thus, once the separations between every pair of peaks sat- 



isfies \qi — qj\ > D, equations ( [L.18| ) will become essentially Pi = and = pi 



Subsequently, the peak positions of the separated pulsons will undergo free 
linear motion pi = Ci, q^ = Cit + g° for a set of 2N constants, the speeds q 
and phases for z = 1, . . . , A^. Thus, the peaks will separate linearly in 
time, in proportion to their difference in heights, and become linearly ordered 
according to height, with the faster ones to the right. 

Numerically, the finite-dimensional invariant manifold of superposed solu- 
tions given in ( 1.15 ), (|1.16 ) and satisfying ( 1.18| ) shall be shown for various 
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confined initial velocity distributions and choices of pulse shape g{x) to de- 
scribe stable pulses which interact elastically like solitons do, sort themselves 
according to height (as expected for geodesic motion), and dominate the solu- 
tion of the initial value problem via their two-body interactions. 



2 Pulson-Pulson and Pulson- Antipulson Inter- 
actions 

Before we show numerical solutions of initial value problems for various choices 
of g{x), we shall present the solution of the 2-pulson and pulson-antipulson 
interaction equations for arbitrary g{x). We follow [||, for the analytical 
solution of the two-peakon and peakon-antipeakon interaction dynamics for 
the Camassa-Holm equation. For N = 2, the collective Hamiltonian ( |1.17| ) 
becomes 

H =^{pl+pl)+PiP2g{q2-qi). (2.1) 

Defining sum and difference canonical variables as P = pi -|- Q = Qi + 12, 
and p = P2 — Pi, Q = 12 — Qi, respectively, transforms this Hamiltonian to 

H = lp'-\{P'-p'){l-g{q)), (2.2) 

which is independent of the sum coordinate Q, so that P is conserved. Thus H 



and P are both constants of motion. Hence, equation ( p.2| ) relates the phase 
space coordinates {p, q) by 

, 4:H-2P^ . -4ClC2 . ^o . s 

upon writing the values of the constants H and P as H = |(cf + C2) and 
P = (ci + C2), for asymptotic speeds Ci and C2. In this relation, the condition 
g{0) = 1 for g = will produce singular behavior for p, should the pulsons 
overlap. For 2-pulson, or 2-antipulson coUisons, we have C1C2 > 0, so the peak 
separation q cannot vanish in these cases for real p. However, for pulson- 
antipulson coUisons, we have C1C2 < 0; so g = may occur and the relative 
momentum p diverges when this happens. 
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The sum and difference variables obey the canonical equations 

P=-2|| = 0, Q = 2% = P{l + g{q)), 

p=-2^ = -UP' - V')g\q) , g = 2|f = p(l - g{q)) . 



(2.4) 



dq 2V" ^ 'if ' 'i " dp 

Conservation of the total momentum P = YliLiPi holds for arbitrary A^, and 
for the case = 2 it is sufficient for solvability of the dynamics. Substituting 
p from ( p.3| ) into the q equation in ( p.4[ ) gives a quadrature formula for the 
dynamics of the separation between peaks, q(t). Namely, 

±t-to = / YFi- 2.5 

A(io) [PV(g') - 4i^^?(gO + (4^ - ^ 

Seen as a collision between two initially well-separated "particles" with initial 
speeds c\ and C2, the separation q{t) reaches a nonzero distance of closest 
approach for 2-pulson collisions with c\C2 > and this quadrature formula 
produces only a phase shift (a delay or advance of position relative to the 
extension of the incoming trajectory, without a change in asymptotic speed) of 
the two asymptotically linear q{t) trajectories. These trajectories are shown for 
several choices of g{x) in Figure pTT|. The phase space trajectories of rear-end 



collisions for three typical Green's functions are shown in Figure |2]^. Notice 
that p is nonsingular and q remains positive throughout, so the particles retain 
their order. The momentum transfer in these rear-end collisions occurs rather 
suddenly over a small range of separation distance q, especially for wave forms 
with compact support. 

Head-on pulson-antipulson collision. For the special case of completely 
antisymmetric pulson-antipulson collisions, for which pi = —p2 = —p/2 and 
qi = —q2 = —q/2 (so that P = and Q = 0) the quadrature formula ( |2.5| ) 
reduces to 

1 r^*) dq' , , 

±{t-to) = ^= , ^ 2.6 

For this case, the conserved Hamiltonian (|2.2|) is 



H = ^[l-giq)]. (2.7) 
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10 15 20 25 30 35 



qi,q2 

Figure 2.1: Space-time trajectories qi(t) and q2{t) for three typical Green's 
functions. For any spatially confined Green's function, these space-time plots 
become linear, asymptotically in time. Note that the relative separation q = 
q2 — qi remains positive, so the particles retain their order. 



After the collision, the pulson and antipulson separate and travel oppositely 
apart; so that asymptotically in time 5f(g) — > 0, p — >• 2c, and H ^ c^, where c 
(or — c) is the asymptotic speed (and amplitude) of the pulson (or antipulson). 
Setting H = in equation ( |2.7| ) gives a relation for the pulson-antipulson 
{p, q) phase trajectories for any Green's function, 

2c , , 

P = T^===. (2.8) 

Notice that 1/p passes through zero when g 0, since g{Q) = 1. The relative 
momentum p, initially at —2c, diverges to — oo at the "bounce" collision point 
as g ^ 0+. Then p changes sign, q increases again and p asymptotes to 2c from 
above. Thus, p diverges as g ^ 0^ and switches branches of the square root, 
from negative to positive. Note also that g > throughout, so the particles 
retain their order. The phase space trajectories of head-on collisions for three 
typical Green's functions are shown in Figure p.3| . 

In terms of relative momentum and separation, p and g, the solution ( |1.16| ) 
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0.6 




Figure 2.2: Phase space trajectories p{t) vs q{t) of rear-end collisons for three 
typical Green's functions. For rear-end colhsons these trajectories are nonsin- 
gular. Note that p = P2 — Pi changes sign, while q = q2 — Qi remains positive; 
so the particles exchange their momentum, but retain their order. 



for the velocity u in the head-on pulson-antipulson collision becomes 



— P P 
t) = -i^9{x + q/2) + -g{x - q/2) . 



(2.9) 



Using equation ( |2.8| ) to eliminate p allows this solution to be written as a 
function of position x and the separation between the pulses q for any pulse 
shape g{x) as 



u{x,q) 



g{x + q/2) - g{x - q/2) 



:2.10) 



where c is the pulson speed at sufficiently large separation and the dynamics 
of the separation q{t) is given implicitly by equation ( p. 61 ) with \/ AH = 2c. 
Equation (|2.10|) is the exact analytical solution for the pulson-antipulson col- 
lision. 
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g(q) = e-'^^ 
g(q) = 1-|q| 



"0 1 2 3 

q 

Figure 2.3: Phase space trajectories of head-on coUisions for three typical 
Green's functions. As g — > 0"^ at the coUison point, 1/p passes through zero. 
Note that q remains positive throughout the colhsion, so the particles retain 
their order. In particular, they do not "pass through" each other. 

3 Numerical Simulations 

We shall summarize our choices of Green's functions and briefly discuss the nu- 
merical method used to compare their pde interaction dynamics in numerical 
simulations. These numerical simulations demonstrate the qualitative similar- 
ity in solution behavior for these Green's functions in both rear-end (2-pulson) 
and head-on (pulson-antipulson) collisions. The numerical results show that 
the dynamics of the initial value problem is dominated by Green's function 
wave forms that collide elastically with each other and sort themselves accord- 
ing to height, as solitons do. The reversibility of the numerical dynamics for 
the initial value problem is also shown, by demonstrating that the dynamics 
of a set of separated peaked solitary waves may be reversed to reconstruct a 
smooth initial distribution. 
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3.1 Choice of Green's Functions 

The choice g{x) = e"'"^' satisfies (1 — d'^)g{x) = 25{x), for which the Hamilto- 



nian 



h3i ) becomes 

Hi = - J mu dx = - J {u^ + u1)dx . (3.1^ 



This is the Hi norm which generates the peakon dynamics studied in [§], 0, ||T^, |1T8| 
It is natural to consider higher- derivative norms of this type, such as 

H, = \ f {u^ + ul)dx, (3.2) 



2 „ 

with Us = d^u/dx'^. The corresponding normahzed Green's function 



^E-.e-^^N, (3.3) 



g^'\x) 

with Tj = exp (j — ^))) j = 0, s — 1 and R = ^^Ig^j, satisfies 

(l + (-9^)^)^(^)(x) = |5(x). (3.4) 

The Green's function g^{x) has a discontinuity in its s-th derivative. The pulse 
shape (|3.3| ) is plotted for s = 1 peakons and s = 2, 3 pulsons in Figure [3^] . 
In the high wave number limit, the Hamiltonian ( p.2| ) becomes 

Hs,k^oo = 7: I uldx, (3.5) 



2 „ 

whose corresponding Green's function satisfies 

i-dygi'\x) = 2sS{x). (3.6) 



The Green's functions corresponding to ( p.6|) having compact support are given 
by 

.is)^^, _ / 1 + ^x'^~hgn{x) if |a:| < ( I^^^^ 



g(f){x) ={ (2^^ "^"^^^ ^ y^^J ' (3.7) 

otherwise. 
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Figure 3.1: Pulson shapes for s=l,2, and 3. 




-0.2 
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Figure 3.2: Compacton shapes for s=l,2, and 3. These shapes widen as s 



increases. 



These are depicted in Figure according to ( |3.7D for s = 1,3, but the opposite 
convexity is chosen for s = 2. These "compacton" Green's functions serve as 
examples of travehng wave solutions to (|1.1| ) with compact support. 

In addition to the aformentioned Green's functions, we also test a Gaussian 
and a composite Green's function consisting of one central s = 2 compacton 
and two outlying s = 1 compactons. 
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3.2 Numerical Method 



Equation (|1.1| ) is advanced in time with a fourth order Runge Kutta scheme 
whose time step is chosen (by numerical trial and error) so as to ensure less 
than a 1% drop in amplitude over 8 periodic domain traversals. We keep track 
of the evolution of the Fourier modes of m{x, t) and compute the spatial deriva- 
tives in ( |1.1| ) pseudospectrally. To transform between u{x,t) and m{x,t) we 
convolve m{x,t) with the Green's function in Fourier space as m„ = Ngnrhn, 
where N is the number of Fourier modes, which is kept at 2048 for each run. 
This method of transforming between m and u proves convenient for arbi- 
trary Green's functions because the transformation can be computed numer- 
ically and the relationship between m and u does not need to be determined 
explicitly. Because antisymmetric perturbations to the zero solution are un- 
stable [Q, 0, and the numerical approximation of the nonlinear terms have 
aliasing errors in the high wave numbers, we apply the following high pass 



filtered artificial viscosity pi 



if \k\ < f 

m = { 2em - 1) if ^ < < , (3. 



N ^1 " 16 
£ if \k\ > 



3N 



where e = 0.01 for the present simulations. 



3.3 Initial Value Results 

Beginning with a normalized Gaussian distribution u{x, 0) = ^r^e^*^^'^^"^^/""^ 
we first show the behavior of the initial value problem for the integrable cases 
of peakons and s = 1 compactons. From the asymptotic speeds of 

the emergent integrable pulses can be determined as the eigenvalues of the 
Sturm-Liouville problem 



CH 



2A 



m{x, t) 



^ = 0, 



(3.9) 



where m{x, t) is the momentum density associated with the velocity of the 
initial distribution and 5ch is 1 for the peakons and for the s = 1 compactons. 
Camassa et al. solve (|3.9|) analytically with m(x, 0) = Asech^{x). The 



relation of equation ( |3.9| ) to the spectral problem for the Schrodinger equation 
in the case of positive "potential," m(x, t) has recently been discussed in Beales 
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et al. |25]. See also Constantin |^| for a similar discussion of the spectral 



problem for equation ( p.9|) in the periodic case. 

In the simulations we investigate here the asymptotic speeds are given by 
the heights of the emergent pulses computed numerically from ( p. . 1|) . Figures 



T7\\ and [7.2| depict the emergence of peakons and s = 1 compactons from 
the initial Gaussian distribution. The integrable behavior is evidenced as the 
pulses in both cases collide elastically as they recross the periodic domain. 



Note the presence of roughly nine distinct s = 1 compactons in Figure 7.2 



versus three peakons in Figure |7.1| . While the maximum speed is identical for 
both cases, the number of emergent pulses increases for pulses of smaller area. 

In Figure TA we see the interaction of the larger peakons with the slower 
residual distribution that is left behind by the initial distribution after some 
early evolution has occurred. Notice that the last two peakons emerging in 
the upper right hand corner of the contour plot experience a positive phase 
shift (to the right) shortly after the sixth peakon emerges. This is consistent 
with the increase in amplitude of the peakon as it grows out of the initial 
distribution and propagates away, enhanced by the interaction of the larger, 
faster peakons with the much slower residual distribution, whose peakon con- 
tent has not yet emerged. This rather complex interaction in which the sixth 
peakon experiences a significant amount of growth as it traverses the upper 
right hand corner of the figure could have been avoided by removing the resid- 
ual distribution remaining after the first few peakons had emerged. However, 
the Gaussian initial distribution for this figure is comprised of a countable 
infinity of peakons interacting nonlinearly among themselves. Thus, removing 
the residual distribution this way would alter the intended initial value prob- 
lem perhaps unacceptably. Should we desire only a certain finite number of 
peakons, say, six of them, to emerge from an initial distribution, then we could 
use the inverse problem to determine what initial distribution is formed from 
this number of pulses. We discuss the inverse problem as well as the reversible 
nature of equation ( p..l| ) later in the paper. 



Next, we discuss the nonintegrable results in Figures 7.3, 7.4, 7.5 and 



[7.6| . In Figure we see s = 2 pulsons emerging from an initial Gaussian 
velocity distribution. Note the curvature of the space-time trajectories in the 
contour plot for the third and fourth pulsons. This occurs because the larger 
width of the s = 2 pulsons compared to the peakons in Figure |7.1| implies a 
longer interaction period between pulson pairs. Again we see the effects of the 
interaction of the larger pulsons with the residual distribution as they recross 
the periodic domain. Most notably we see the continued emergence of pulsons 
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on the "upwind" side of the initial distribution, leading to the emergence of 



a pulson at the top center of the contour plot in Figure [7^^. Figures and 
|7.5| also depict this phenomenon of curvature in the space-time trajectories 
during the emergence of the s = 2 compactons and Gaussons. In Figure 
|7.4| the curvature in the space-time trajectories is not only because of the 
prolonged interaction between the s = 2 compactons, but also because of their 
interaction with the residual initial distribution as they recross the domain. 
We do not show the Gaussons emerging from the initial s = 1 compacton 



shape and colliding with the initial residual distribution in Figure ^]6| because 
there are too many pulses to make a clear plot. However, this figure helps to 
support our main conclusion: that an arbitrary (even) wave form given by the 
Green's function g{x) can be made to emerge from an (essentially) arbitrary 
confined initial condition. In this case we have chosen a wide s = 1 compacton 
for the initial distribution partly because its shallow discontinuity imposes 
less stringent conditions on the numerics, and also to demonstrate the effects 
that the width of the initial condition has on the number of emergent pulses. 
Namely, the width of the initial distribution must be greater than the single 
pulson width prescribed by the function g{x), and wider initial distributions 
produce a larger number of pulsons. 

In summary. Figures fflSj [7^ , |775| and |7]^ demonstrate that: (1) the shape of 
each wave form emerging in the initial value problem is controlled by the choice 
of g{x)\ and (2) the dynamical behavior of the emergent pulses is qualitatively 
the same as the integrable behavior of the peakons and s = 1 compactons. 

3.4 Rear-End Collisions 

We compute the evolution governed by equation (|1.1|) with two wave forms (for 
a particular choice of Green's function) located initially at x = and x = L/2 
with amplitudes ci = 1 and C2 = 1/2, respectively. This situation is arranged 
to investigate the nature of rear-end collisions both for the integrable and 
nonintegrable cases. Beginning with the integrable cases depicted in Figures 
|7[^and|7]^, we see that, consistent with the findings of the slower peakon in 



Figure 7.7 experiences no phase shift, but the faster peakon experiences a shift 



to the right. This occurs after each collision, thus causing the second collision 



to occur to the left of the first. In Figure |7.8| we see that the faster s = 1 
compacton trajectory experiences a phase shift to the right while the slower 
one experiences a shift to the left. Thus, perhaps not unexpectedly, the phase 
shift of the interaction depends on the pulse shape. The sudden leftward phase 
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shift of the slower trajectory indicates that momentum is being transferred 
rightward rather quickly in the collision process across the entire width of 
both wave forms. This is something like the sudden transfer of momentum 
in the collision of billiard balls. For this reason, the magnitude of the phase 
shift increases with pulse width. Note that the pulsons exchange momentum 
in each collision, rather than passing through each other. 

The nonintegrable cases depicted in Figures fO| , [7.10| , [7.11| , and [7.12| all 



show similar elastic collisions resulting only in phase shifts. In each case the 
faster pulse trajectory experiences no shift, or a shift to the right, while the 
slower one experiences a shift to the left. The magnitude of the shift increases 
with pulse width. The pulse width is the largest for the multicompactons in 
Figure [7.12| and they experience the largest phase shift of the four noninte- 



grable cases shown here. The Gaussons in Figure |7.11| , on the other hand, 
experience the least phase shift and they are the narrowest of the four nonin- 
tegrable pulsons we treat. 

3.5 Head-on Collisions 

We use the exact solution u{x, q) for arbitrary Green's functions in equa- 
tion ( p.lO|) to determine the velocity profiles in head-on antisymmetric pulson- 



antipulson collisions. The complexity of the wave shapes and the strengths of 
the various singularities which form according to equation ( |2.10|) during these 
interactions for many choices of pulson shapes are beyond the capabilities of 
most numerics. For example, the spatial derivative of equation (|2.10|) at x = 
gives 

u,{Q,q) = ^JL==g\q/2). (3.10) 

Hence, whether a vertical slope in u forms at x = depends on the choice of 
the Green's function g{x). Colliding peakon-antipeakon profiles are depicted 
in Figure [7.13| . As discussed in the solution u in this case tends to zero as 



g ^ 0, as the slope limg^o "^xlO, g) diverges to minus infinity. Also, p — oo 
as g — 0, so as to maintain constant energy and zero total momentum in the 
particle system. This divergence phenomenon occurs similarly for the s = 1 



compactons, as shown in Figure [7!T^ . 



As with the peakons, the colliding s = 1 compactons develop a vertical 
slope in velocity at a; = when q = 0. Upon making initial contact when 
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g = 2, the two triangular pulses develop a maximum negative slope at a; = 
of "^^(0, g)|gg[i^2] = ~2c which remains constant until g = 1, whereupon 
the two pulses begin to clip their peaks and become trapezoidal as the slope 
at a; = tends toward minus infinity via limg^o+ ^x(0, g) = —2c/ \/\q\, cf. 
equation (|3.1CI| ). When g = 0, the two colliding pulses diminish to m = and 



then "bounce" apart to reverse the aformentioned process, so that M^(0,g) = 
+2c/ ^/\q\ until g = 1 with Mx(0, g)|gg[i,2] = +2c. Thereafter, the reformed 
pulses separate in opposite directions. Note that g approaches zero, but does 
not change sign; so the "particles" with phase space coordinates pj, g^, i = 1,2, 
keep their order, just as for the rear-end collisions. 

The most significant difference between the integrable and nonintegrable 
pulson-antipulson collisions is that the solution u{x, g) does not necessarily 
tend toward zero as g ^ for the nonintegrable cases. Another difference is 
that a verticality in slope does not necessarily develop at x = 0. Moreover, 
such a verticality may not develop at all. As an example of a nonintegrable 
pulse that does not approach -u = as g — 0, but does develop two verticalities 
at that instant, we consider a parabolic pulse such that 

(X) = I < 1, 

^ I otherwise. 

The antisymmetrically colliding parabolas are shown in Figure |7.15| . As the 
parabolic pulson collides with its antipulson, two verticalities develop at x = ± 1 
rather than at x = 0. By equation (|2.1CI|) , this interaction creates a limiting 
distribution such that 

lim^M(x, g) <^ ^ otherwise. (^•-'-2) 



As an example of a nonintegrable pulson that does not tend toward u = 
and develops no verticalities in slope as g ^ 0, we consider the antisymmetri- 
cally colliding Gaussons depicted in Figure 7.16 . Here the limiting distribution 
is 

lim u(x,q) = -2cxe-''\ (3.13) 

When an antisymmetric collision results in a limiting distribution as in equa- 
tions ( |3.12| ) and ( p.l3|) as g — 0^, the solutions must flip instantaneously as 



they approach g = 0^, and bounce apart so that the post-collision interaction 
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will be identical but opposite in polarity to the pre-collision interaction. This 
situation also occurs similarly for the antisymmetric collision of s = 2 pulsons, 
shown in Figure [7.1 7| . However, the analytical solution for the limiting g ^ O"*" 



distribution is not so easily expressed in this case, so we omit it here. 

An even more complex situation develops when two s = 2 compactons 
collide antisymmetrically, as shown in Figure [7!T8| . In this case, not only do the 
colliding pulses develop two verticalities, but there is no limiting distribution 
as g — s> 0~^. In fact, one can show for s = 2 compactons (with the factor of 1/3 
omitted for convenience) that 

u{x = ±1, q) = fl + (3.14) 



so that the velocity at x = ±1 becomes unbounded as g 0+. However, 
as soon as q reaches zero, the unbounded solution for u switches sign and 
the pulses restore themselves with reversed polarity and move apart. As a 
final example we show in Figure |7.19| the antisymmetric collision of "multi- 



compactons," again with the factor of 1/3 omitted from the central hump for 
convenience. Here we use a "waterfall" plot to depict the collision since an 
overlay would be too confusing. This collision results in the creation of eight 
verticalities at x = ±1, ±2, ±3, and ±4. Again, as g ^ O"*", the solution flips 
polarity instantaneously and the interaction reverses as the multicompactons 
move apart. 



3.6 Reversibility 

The results in Section |3.3| for the initial value problem show that a train of 
pulses in the prescribed pulse shape g{x) emerges from any initial distribution 
and dominates the initial value problem. In this process, the initial distribu- 
tion breaks up into a discrete number of pulsons, each of which travels at a 
speed equal to its height, and the pulsons colhde elastically among each other 
as they recross the periodic domain. An interesting feature of equation (|1 . 1| ) is 
that it is time reversible, that is, invariant under t — > — t and u —u. Thus, 
replacing dt with —dt at any point in the numerical simulations causes the pul- 
sons emerging in the initial value problems to reverse their evolution sequence 
(but not their heights) and collapse together into the initial distribution. Thus 
the time-reversed series of collisions reforms the original initial distribution at 
t = 0. This process enables us in principle to determine the parameters pj. 
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gi, i = 1, . . . ,N, on the invariant manifold by running the pde ( |1 . 1| ) forward 
in time until pulsons have emerged, then reversing the solution to recre- 
ate the original smooth initial distribution. Similarly, we may run the pde 
forward until time T, determine the parameters {pi, qi) at this time, remove 
any residual velocity distribution from which pulsons have not yet emerged, 
then time-reverse the ordinary differential equation (ode) dynamics for a pe- 
riod T to create an initial distribution containing only the pulson content that 
will emerge before time T. Should this process, denoted in obvious notation 
as [ODE{—T) o o PDE{T)]u{x^Q), produce an accurate approximation 
un{x, 0) of the original distribution u{x, 0), it would be a useful procedure for 
approximately determining the pulson invariant manifold content of a given 
initial velocity distribution. 

Figures |7.2CI| and [7. 2 1| were constructed by running the pde forward in time 
from an initial Gaussian distribution until trains of three peakons in Figure 
7.20\ and nine s = 1 compactons in Figure [7.21| emerged. At this point time 
was reversed, and the pde was run in reverse for twice as long, then the 
solution was reflected in space. In this process trains of pulsons in reversed 
order (with the tallest ones behind) reassembled moving rightward into the 
original initial condition and then re-emerged as before (with the tallest ones 
ahead). The same "geodesic pulson scattering" figures could be generated by 
running the initial Gaussian distribution both forward and backward in time, 
then superposing the results. Formally we have S{t)RxS(t) = Id (where S(t) 
is pde evolution for time t and is reflection in space), for initial conditions 
that are symmetric about the origin. Thus these figures could also have been 
generated by the process of evolution, reflection, then evolution again. This 
time reversal symmetry property indicates that the initial distribution may be 
regarded not just as a sum of pulsons, but as a moment in time at which a set 
of pulsons moving geodesically has collided to form a smooth distribution, in 
this Gaussian. 



4 Conclusions 



In equations (|1.1| ) and (|1.2|) we have introduced a new family of hyperbolic 
pde describing a certain type of geodesic motion and for some members of this 
family we have studied its particle-like solutions. We call these solutions "pul- 
sons." We have shown that, once initialized on their invariant manifold (which 
may be finite-dimensional), the pulsons undergo geodesic dynamics in terms of 
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canonical Hamiltonian phase space variables. This dynamics effectively repro- 
duces the classical soliton behavior. We conjecture that this behavior occurs 
because of the preponderance of two-body elastic-collision interactions for the 
situation we consider of confined pulses and confined initial conditions. 

The dynamics we study in this framework of geodesic motion on a finite- 
dimensional invariant manifold seems to account for all of the classical soliton 
phenomena, including elastic scattering, dominance of the initial value problem 
by confined pulses and asymptotic sorting according to height - all without 
requiring complete integrability. Thus, complete integrability is apparently 
not necessary for exact soliton behavior. Regarding the "formation" of the 
pulsons: in fact, the pulsons must always be present, since they compose an 
invariant manifold. One discerns them in a given (confined) initial condition 
by using the isospectral property in the integrable cases, or just by waiting for 
them to emerge under the pde dynamics in the nonintegrable cases. However, 
there is no "pattern formation" process in this dynamics. Rather, there is an 
"emergence of the pattern," which is the pulson train. 

This conclusion is illustrated by Figures [7.20| and |7.21| showing geodesic 
scattering of an incoming set of pulsons that collapses into a smooth "initial" 
distribution, then fans out again into an outgoing set of pulsons that is the 
mirror image of the incoming set. This parity-reflection is an implication of 
reversibility, as well; since reversing the dynamics of a set of separated pulsons 
reassembles them into a smooth "initial" distribution, which is just their sum, 
with appropriate values for their initial "moduli" parameters pj(0),gj(0), i = 
1, . . . ,N. The moduli parameters pi, qi are collective phase space coordinates 
on an invariant manifold for the pde motion. Once initialized, these collective 
degrees of freedom persist and emerge as a train of pulses, arranged in order of 
their heights. Moreover, since their velocity is their peak height, small residual 
errors departing from the pulson superposition in the assignment of initial 
parameter values Pi(0), gi(0), i = 1, . . . , A^, do not propagate significantly. On 
the real line, such residual errors are simply left behind by the larger pulsons 
traveling more quickly. And on the periodic interval, these errors provide an 
occasional phase shift perturbation to the stable pulsons. These perturbations 
are seen in the figures as regions of curvature in the space-time trajectories of 
the pulsons. 

Our main conclusions regarding the elastic collisions of the pulsons and 
their role in the pde initial value problem are the following: 



1. Two-pulson interactions are elastic collisions that conserve total kinetic 



Geodesic Soliton Behavior 



23 



energy and transfer momentum, as measured by the asymptotic speeds 
and peak heights of the pulsons. The motion of the pulsons is governed 
by geodesic kinematics of particles with no internal degrees of freedom; 
but with an interaction range, or spatial extent, which is given by their 
Green's function, or wave form, g{x). Thus, these interactions are analo- 
gous to the collisions of billiard balls. This is especially clear for pulsons 
of compact support. 

2. The exact solution in equation ( |2.10|) determines analytically the com- 
plex wave shapes and the strengths of the various singularities that may 
form during antisymmetric head-on collisions of pulsons and antipulsons. 
These singularities are verticalities in slope at which the velocity may 
diverge, take a finite value, or even vanish, depending on the choice of 
Green's function. At the moment of collision, the pulson and antipul- 
son "bounce" apart and the polarity of the antisymmetric wave form 
reverses, u — > — m. Again this is analogous to billiards, but with the 
difference that the reversal of momentum reverses the polarity of the 
pulson velocity profile. Note that the wave forms reverse polarity and 
bounce apart; they do not actually "pass through" each other, although 
they may appear to do so. 

3. The dynamics on the pulson invariant manifold is dominated by the pre- 
ponderance of two-body coUisons. To the extent that the initial value 
problem for the pde takes place on the pulson invariant manifold, these 
two-pulson collisions should also dominate the solution of the pde dy- 
namics. We conjecture this is so for most choices of confined pulse and 
initial conditions of finite extent for the family of pde (|1 . 1| ) with u given 
in ([13). 

Thus, geodesy governs the motion and explains the elastic-scattering soli- 
ton phenomena we observe, without providing (or even requiring there exists) 
a means of analytically solving the initial value problem for this family of pde. 
Even the motion on the pulson manifold is not integrable for > 2, except 
in the cases corresponding to the Camassa-Holm, and Dym integrable pde. 
Similar geodesic but nonintegrable families of equations showing soliton be- 
havior on a finite-dimensional invariant manifold may exist for other integrable 
geodesic equations, such as KdV. 
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6 Appendix on Numerical Stability 



As discussed in Section |3.2| , we add dissipation in the higher modes for the 
simulations covered in this paper because the Fourier pseudospectral method 
used is unstable to perturbations at the grid scale. These perturbations arise 
as a result of the Gibbs phenomenon which causes an overshoot to occur at 
points of discontinuous slope in u. Such an overshoot forms negative values 
of the function m, which can artificially create antipulsons that emerge and 
interact with positive pulsons to form verticalities in slope as described in 



Section |3.5| . As an example, we tested the stability of the numerical method 



by initializing the domain with a perturbation at a single grid point, as 

r 10-3 lit = N/2, 
' 1^ otherwise, 

for i = 1, . . . ,N . This initial distribution was unstable and caused exponential 
growth in the neighboring grid points, such that Un/2-i — * +oo and Um/2+i 
— oo, resulting from the creation of a pulson-antipulson pair whose peaks grew 
unstably. We then ran the same initial condition with a dissipation coefficient 
oi e = 0.01 in equation ( |3.8| ), which stabilized it. We also ran it again with 
a nonzero linear dispersion coefficient oi k = 0.01 in equation (|1 . 1 1| ) which 
stabilized it as well. Adding dissipation damps the higher modes as specified 
in Section |3.2| while adding linear dispersion smooths discontinuities in m, thus 
effectively damping the higher modes. Thus, either dissipation or dispersion 
will suppress instability due to perturbations at the grid scale. 

Figure ^]T] compares dissipative and dispersive dynamics for the case of 
s = 1 compactons emerging from an initial Gaussian velocity distribution. 
Linear dispersion (k ^ 0) acts to smooth out discontinuities as well as hasten- 
ing the emergence time of each pulson. This is demonstrated by the separation 
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Figure 6.1: s = 1 compactons emerging from a Gaussian of width a = 3 
and unit area centered about x = 10 on a periodic domain of length L = 50 
demonstrating the effects of the hnear dispersion coefficient n. The sohd hne 
depicts the k = case with added dissipation while the dashed line depicts 
the K = 10~^ case with no dissipation. 

between the leading pulses for the n = 10~^ and k = cases. Despite the neg- 
ligible difference in speed between the two, the leading n = 10~^ pulson is far 
ahead of the leading k = pulson, because the former was emitted sooner. 
Stability is acquired by suppressing the higher modes, as shown in Figure |6.2i 
Here the power spectrum of one compacton is shown for both values of k. 
Introduction of linear dispersion with k = 10"^ damps the higher modes by 
roughly three orders of magnitude. 




-100 -50 50 100 

k 



Figure 6.2: Power spectra of one pulse for the k = and n = 10 ^ cases. 
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Figure 7.1: s = 1 peakons emerging from a Gaussian of unit area and a = 5 
centered about x = 20 on a periodic domain of length L — 100. A secondary 
peakon emerges in the upper right hand corner of the contour plot as the faster 
peakons recross the domain and interact with the residual distribution. 
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Figure 7.2: s = 1 compactons emerging from a Gaussian of unit area and 
(7 = 5 centered about x = 20 on a periodic domain of length L = 100. More 
compactons emerge here than peakons in Figure because the compactons 
have smaller area. The residual distribution around x = 20 causes the space- 
time trajectories to curve as the larger, faster s = 1 compactons interact with 
the slower ones that comprise the residual. 
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Figure 7.3: s — 2 pulsons emerging from a Gaussian of unit area and a — 5 
centered about a; = 20 on a periodic domain of length L — 100. The space- 
time trajectories curve just after the s = 2 pulsons emerge from the initial 
distribution because their large width enhances their interaction range. Note 
the phase shift in the space-time plot occurring at the trailing edge of the 
residual distribution as the larger, faster pulsons recross the domain. 
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Figure 7.4: s — 2 compactons emerging from a Gaussian of unit area and a — 5 
centered about a; = 20 on a periodic domain of length L = 100. The space- 
time trajectories curve shghtly as they emerge from the initial distribution 
and then straighten out when the compactons are no longer interacting. The 
space-time trajectories curve again as the compactons recross the domain and 
interact with the residual distribution. 



Geodesic Soliton Behavior 



33 




Figure 7.5: Gaussons emerging from an initial Gaussian of unit area and a — 5 
centered about x = 20 on a periodic domain of length L = 100. The interaction 
of the faster Gaussons with the residual "ramp" distribution causes another 
Gausson to emerge at the top center of the contour plot. The interaction 
between nearby Gaussons and between Gaussons and the residual distribution 
causes their space-time trajectories to curve. 



Geodesic Soliton Behavior 



34 




10 20 30 40 50 60 70 80 90 100 



X 

Figure 7.6: Gaussons emerging from an initial s = 1 compacton of unit area 
and width 40 centered about a: = 20 on a periodic domain of length L — 100. 
In an effort to avoid confusion, the Gaussons are not shown recrossing the 
domain. 
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Figure 7.7: Rear-end collision dynamics for s = 1 peakons. The faster peakon 
moves at twice the speed of the slower one. For this case, both coUisions result 
in a phase shift to the right for the faster space-time trajectory, but no phase 
shift for the slower one. The phase shift causes the second collision to occur 
slightly to the left of the first. 
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Figure 7.8: Rear-end collision dynamics for s = 1 compactons. The faster 
space-time trajectory experiences a phase shift to the right while the slower 
one experiences a phase shift to the left. 
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Figure 7.9: Rcar-cnd collision dynamics for s = 2 pulsons. The faster space- 
time trajectory experiences a phase shift to the right while the slower one 
experiences a phase shift to the left. The size of the interaction region is 
proportional to the sum of the pulse widths. 
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Figure 7.10: Rear-end collision dynamics for s = 2 compactons. The faster 
space-time trajectory experiences a phase shift to the right while the slower 
one experiences a phase shift to the left. 
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Figure 7.11: Rear-end collision dynamics for Gaussons. The faster space- 
time trajectory experiences a phase shift to the right while the slower one 
experiences a phase shift to the left. Because they are so narrow, Gaussons 
have the smallest phase shift of the four nonintegrable cases depicted in this 
figure and Figures [^9] , [7.10| , and |7.12 . 
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Figure 7.12: Rear-end collision dynamics for multicompactons. The faster 
space-time trajectory experiences a phase shift to the right while the slower 
one experiences a phase shift to the left. The large width and hence large 
interaction range of the multicompactons causes them to have the largest phase 
shift of the four nonintegrable cases depicted in this figure and Figures [7.9| , 
7l0| , and 
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Figure 7.13: Head-on collision dynamics for s = 1 peakons at five pre-collision 
separations q. The velocity tends to zero as the separation between the peakons 
tends to zero and a verticality in slope develops at a; = 0, as g — > 0+. At 
this collision point, the solution reverses its polarity by flipping across the 
horizontal axis (so that u — > —u) and the peakons move apart in opposite 
directions. 
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Figure 7.14: Head-on collision dynamics for s = 1 compactons at five pre- 
collision separations q. As the compactons meet, the slope at a; = remains 
at Ux{0,q) = —2 until their peaks are clipped when q = 1. Thereafter the 
slope at a; = begins to increase and eventually becomes vertical when g = 
just as the velocity vanishes. At this moment the solution reverses polarity 
and the compactons re-emerge, moving apart in opposite directions. 
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Figure 7.15: Head-on collision dynamics for parabolas at five pre-collision 
separations q. Parabolas collide to form a finite distribution when g = 0+ and 
two verticalities in slope appear at x = ±1. After the collision the solution 
reverses polarity and the parabolas move apart, preserving the antisymmetry. 
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Figure 7.16: Head-on collision dynamics for Gaussons at five pre-collision sep- 
arations q. Gaussons collide to form a finite distribution when g = 0+ and 
develop no verticalities in slope. After the collision the solution reverses po- 
larity and the Gaussons move away from each other. 
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Figure 7.17: Head-on collision dynamics for s = 2 pulsons at five pre-collision 
separations q. As for the Gaussons in Figure [7.16| , s = 2 pulsons collide to 
form a finite distribution when g = O"*" and develop no verticalities in slope. 
At the moment of collision, the solution reverses polarity and then the s = 2 
pulsons move apart antisymmetrically. 




Figure 7.18: Head-on collision dynamics for s = 2 compactons at five pre- 
collision separations q. The s = 2 compactons collide to form two verticalities 
in slope at a; = ±1 as well as unbounded growth in u{x, q)\x=±i as q ^ 0+. At 
this moment, the solution reverses polarity and the compactons move apart in 
opposite directions. 
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Figure 7.19: Head-on collision dynamics of multicompactons depicted at sep- 
arations of g = 10, 5, 4, 3, 2, 1, and 0.1. We avoid confusion and show a 
"waterfall" plot to depict the multicompactons rather than an overlay. Here, 
verticalities as well as unbounded growth in u{x, q) develop at x — ±1, ±2, 
±3, and ±4 as — > 0+. The last distribution is scaled by a factor of 1/2 to 
keep the large amplitudes resulting from the approaching unboundedness of 
u{x,q) within the plot. 




Figure 7.20: s = 1 peakons colliding to form a Gaussian and then reemerging 
symmetrically thereafter. 
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Figure 7.21: s = 1 compactons colliding to form a Gaussian and then reemerg- 
ing symmetrically thereafter. 



